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C/3 \ ABSTRACT 

G\ ' 

The variability in the infrared to millimetre emission from microquasar GRS 
. 1915+105 is believed to be dominated by the system's relativistic jet. In this paper 

J> ■ we develop a time-dependent version of the Blandford & Konigl (1979) jet emission 

t-H | model and apply it to the oscillations in the infrared and millimetre emission from 

! GRS 1915+105 observed by Fender & Pooley (2000). The resulting model provides 

' a reasonable description of the observed flux oscillations from GRS 1915+105. From 

a fit of the observed time lag between the flux peaks in the infrared and millimetre 
^— ^ ■ emission together with the flux normalisation we were able to determine the model 

' parameters for the GRS 1915+105 jet. We find that to achieve the observed flux lev- 

els with the model requires an unphysically large electron density within the jet. We 
' therefore conclude that the Blandford & Konigl (1979) model cannot explain these 

Q-f observations, either because it does not provide the correct description of the emission 

q . from microquasar jets, or because the observed emission variations do not originate 

' in the jet. 
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1 INTRODUCTION This suggests a flat spectrum over a frequency range that 

encompasses three orders of magnitude. Although flat spec- 

The microquasar GRS 1915+105 snows time variability m , , . , , , , r ^rr, , inr ,i 

tra had previously been observed from GRa 1915+105, these 

its emission at all observable wavelengths. Its high energy ' e , , 

, measurements were the first to suggest they extend to the 

emission (X-rays) is believed to be created in the inner radii . „ , 

. , , . ,. , . / t • \ infrared, 
of the system s accretion disc, and its fow energy (radio) 

emission is primarily from synchrotron emitting ejecta lo- 
cated at large radii (r > 10 11 m) within the system's jet. Flat spectra from compact radio sources such as ma- 
in this paper we investigate the cause of the variability in croquasars have traditionally been explained in terms of 
the infrared, and millimetre emission from GRS 1915+105 partially self-absorbed synchrotron emission from the inner- 
which is believed to originate in the inner regions of the most regions of a narrow, conical, relativistic jet (originally 
relativistic jet. by Blandford & Konigl 1979, hereafter referred to as BK79, 
Observations ranging from radio to infrared frequen- but also by Hjellming & Johnston 1988, Falcke & Biermann 
cies have revealed very large amplitude, quasi-periodic os- 1995, Falcke 1996, and Falcke & Biermann 1999). This inter- 
calations in the emission from GRS 1915+105 (Pooley & pretation is supported by radio observations of blobs moving 
Fender 1997; Mirabel et al. 1998; Fender & Pooley 1998; at apparent superluminal velocities outwards from a central 
Fender et al. 2002). In this paper we will concentrate on the elongated radio core (Mirabel & Rodriguez 1994; Rodriguez 
oscillations detected at millimetre and infrared frequencies & Mirabel 1999; Fender et al. 1999; Dhawan, Mirabel & 
by Fender & Pooley (2000). These oscillations are charac- Rodriguez 2000). Detection of polarised emission strongly 
terised by the short time lags between flux peaks and a flux suggests the synchrotron nature of the emission (e.g. Fender 
ratio, between these two frequencies, that is close to unity. et al. 2002). The temporal and spectral behaviour of the 

emission coming from the blobs is consistent with expecta- 
tions from internal shocks in a relativistic jet flow (Kaiser, 

* e-mail: rsc@astro.soton.ac.uk Sunyaev & Spruit 2000). These synchrotron-emitting blobs 
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are therefore believed to be ejecta from a central compact, 
relativistic jet. 

However, the BK79 model only represents a steady- 
state solution for emission from relativistic jets. Using the 
basic principles of the BK79 model we develop a time- 
dependent variation of this model in an attempt to explain 
the observed infrared and millimetre flux oscillations from 
GRS f 9f 5+105. Whereas previous time-dependent jet emis- 
sion models have only been concerned with the emission 
from the spherical ejecta blobs (e.g. van der Laan 1966, and 
Hjellming & Johnston 1988), the model presented in this 
paper concentrates solely upon the variability in the emis- 
sion from the inner-most radii of the conical jet. In section|5| 
we demonstrate both how the steady-state BK79 model pro- 
duces a flat synchrotron spectrum and how this is dependent 
upon the physical parameters of the jet. The model is then 
adapted to include time variability in section[3] which allows 
a direct comparison between the observations and the model 
in section [I] Finally we discuss the successes and failures of 
the model in explaining the observations in section 



2 THE BLANDFORD-KONIGL 
STEADY-STATE JET MODEL 

Observations of flat spectrum synchrotron radiation are ex- 
plained by the BK79 model as the emission from a super- 
sonic, relativistic jet at radii relatively close to its accelera- 
tion region. The jet has a conical geometry, with a constant 
opening angle, as a consequence of the free adiabatic ex- 
pansion of the jet plasma which is modelled as flowing at 
a constant velocity, Vj. Irregularities in the supersonic flow 
gives rise to shocks that accelerate some of the electrons in 
the plasma to relativistic velocities with a power-law energy 
distribution of 



(1) 



up to a maximum value of 7 max , where 7 is the electron 
energy in terms of the Lorentz factor, k is a normalisation 
value for the number density, and p the power-law index 
(typically ~ 2 from acceleration by shocks, e.g. Bell, 1978). 

Flux conservation requires that the magnetic field 
strength, B, in this conical jet varies as r _1 , where r is the 
radial distance from the origin. The radial dependence of the 
electron density normalisation value is solely determined by 
volume expansion in two dimensions and so k varies as r~ 2 . 
This represents the most important assumption of the BK79 
model for the creation of flat emission spectra. It implic- 
itly proposes the existence of ongoing particle acceleration 
throughout the expanding jet to compensate for the energy 
losses of the relativistic electrons due to adiabatic decom- 
pression. Furthermore, radiative losses due to synchrotron 
emission and inverse Compton scattering of the Cosmic Mi- 
crowave Background (CMB) are also suppressed, allowing 
the maximum electron energy cut-off at 7 max to remain con- 
stant. Thus the power-law distribution of 71(7) holds for all 
times (hence the steady-state jet model) throughout the jet. 
Although, as stated by BK79, this assumption is a require- 
ment for the equipartition of energy between the magnetic 
field and the electrons, the cause of the ongoing particle 
acceleration is unclear as a viable physical mechanism has 
not been identified. Note that this re-acceleration process 



differs from that which is usually attributed to shock-wave 
propagation in jets, in that it must be a continuous, ubiq- 
uitous process that reacts to energy losses in such a way as 
to exactly counter the losses. 

The jet's emission is entirely due to self-absorbed syn- 
chrotron radiation, with emission, e„ (J s" 1 m~ 3 Hz" 1 ST ), 
and absorption, x» ( m_1 ) 5 co-efficients that for a power-law 
distribution of electrons are given by 
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respectively (e.g. Longair 1994). When applied to the jet 
geometry, where k = k(r) and B — B(r), the source function 
varies with radius as 



xAr) 



(4) 



For a jet observed at an inclination of 90° , the path length 
along the line-of-sight is assumed to be equal to the jet's 
half- width, w, which for a conical jet varies as r. Therefore 
the optical depth varies with radius as 



r»(r) = Xr(r)w(r) oc r ' (p+4)/2 v - {p+A)/2 . 



(5) 



To obtain the emission spectrum the equation for ra- 
diative transfer through a homogenous medium is solved, 
which gives a specific intensity of 



Iv — Sv 



(6) 



The observed flux density, F v (W m" 2 Hz" 1 ), is then simply 
calculated as the product of the specific intensity with the 
projected surface area of the jet, 



AF V = I v dfi 
ail = 



(7) 



where D, is the distance to the jet. Hence 



dFv oc r 



3/2 5/2 



-0+0/2 _( p + 4)/2 



dr. 



(8) 



2.1 Analytical solution 

A flat emission spectrum results from integrating the emit- 
ted flux from the base of the jet out to the maximum ob- 
servable radius, r maxo b, where the magnetic field strength 
is too weak to produce synchrotron radiation at a given 
frequency. This radius is dependent upon the critical fre- 
quency, v c (r) oc 7ma X S(r), of those electrons with the max- 
imum Lorentz factor, 7 max , which remains constant as a re- 
sult of the postulated re-acceleration process. If the criti- 
cal frequency of these most energetic electrons falls below 
the observing frequency, v, then the maximum radius from 
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which we receive radiation of frequency v is reached. Hence 
r m „„h oc v~ x for a conical geometry. The flux integration 
may be approximated by separating the optically thin and 
optically thick regions of the jet, 



F„ oc / r v 



3/2 „5/2 v r -(p+4)/2 ^-(p+4)/2 



dr 



+ 



3/2 5/2 i 

r v dr 



(9) 



where the optically thin radius, rth(f) oc u^ 1 , is defined as 
the radius where the optical depth of the jet material for a 
given frequency, v, is equal to unity. This recovers the flat 
spectrum result of BK79 for the steady-state scenario. 

2.2 Numerical study 

When modelling time variability, in J|| the emission region 
is constrained to some fixed length, terminating at a radius 
determined by the properties of the jet flow rather than 
the observing frequency. In this case a flat spectrum can 
also be created, without the requirement that the maximum 
observable radius, r maxo b, lies within the emission region. 
However, the flat region of the spectrum will only extend 
over a finite frequency range, with an optically thick, F v oc 



, slope at the low frequency cut-off, and an optically 



^5/2 

thin, F v oc jy _ ( p_1 >/ 2 ( slope at the high frequency cut-off. 

The flat region of the spectrum encompasses every fre- 
quency for which the jet material, over the length of the 
emission region, is optically thick at the smallest radius and 
has become optically thin at the largest radius. Thus the 
flat region is formed by the summation of self-absorbed syn- 
chrotron spectra created by electron populations at increas- 
ing radii along the jet. As demonstrated in figure Q each 
individual spectrum is shifted towards lower frequencies as 
the electron density and magnetic field strength decrease at 
greater radii. If energy losses due to adiabatic decompression 
are included in the model, then the spectra with lower fre- 
quency peaks (from greater radii) will have decreasing flux 
levels, hence an inverted spectrum is observed rather than 
a flat spectrum. 

Therefore by suitably adjusting the parameters that 
control the optical depth of the jet material as a function 
of radius, the frequency range of the flat region of the spec- 
trum, extending from f-mi-n to v max , may be shifted to agree 
with observations. However, the extent of the frequency 
range, ^max/^mm, is solely determined by the length of the 
emission region. In this section we demonstrate these princi- 
ples by numerically integrating equation |7] which calculates 
the flux over the length of the jet. 

The emission and absorption co-efficients as function 
of radius are calculated from equations |5] and |3] with the 
magnetic field strength, B(r) = Bo(r/ro)~ , and the nor- 
malisation value of the electron density distribution, k(r) = 
ko{r /ro)~ 2 . Here -Bo = B(r — ro), and ko = fc(r = ro). As the 
power-law index of the electron distribution, p, is commonly 
stated to lie between 2.2 and 2.3 for relativistic shock accel- 
eration (e.g. Achterberg et al. 2001), we fix p to a value 2.25. 
The precise value of p is unimportant in practice as when 
varied over the range 2 < p < 3 it only has a negligible effect 
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Figure 1. The formation of a flat spectrum from the summa- 
tion of self-absorbed synchrotron spectra. Here a sample of spec- 
tra emitted from electron populations located at increasing radii 
along the jet reveal the shift of the individual spectra towards 
lower frequencies. The higher the frequency of the spectral peak, 
the closer to the jet origin. 



upon the model results. Hence the source and optical depth 
functions become 



5/2 



J s _1 m~ 2 Hz- 1 sr-\ (10) 



and 



r v {r) = 1.5X 10- 12 fcoB 17/8 ^o (£) ~ 25/S (^) , (11) 

where, for a conical jet, we have w(r) = wo(r/ro). 

The emission spectrum can then be calculated from 
equation |7| by performing a Romberg numerical integration 
(e.g. Press et al. 1992) of 
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(12) 



where all units are in SI unless stated otherwise. For a fixed 
emission region length, and jet opening angle (which deter- 
mines wo), the model has two parameters; ko, and -Bo. To 
simplify the analysis of the model, we introduce a single pa- 
rameter, ao = ko Bq 7 ^ 8 , that uniquely determines the optical 
depth function. If the flux is normalised, then ao becomes 
the only free model parameter, as the other model parame- 
ter, Bo, is determined by the flux normalisation. 

The spectra, normalised to the maximum flux, from an 
emission region of size, r max /r m i n = 10 s , for three different 
values of ao are shown in figure [5] This figure demonstrates 
that an increase in the value of ao, which increases the opti- 
cal depth of the jet material, will shift the spectrum towards 
higher frequencies. An increase in the value of ao will also 
increase the unnormalised flux for all frequencies. 

The low frequency cut-off to the flat region of the spec- 
trum occurs at the frequency at which the jet material at 
r ma x becomes optically thin, so that the entire jet is op- 
tically thick to all lower frequencies. Conversely, the high 
frequency cut-off occurs at the frequency at which the jet 
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Figure 2. Normalised spectra from emission region of length 
fmax = 10 s r m ; n for values ao = 10 2 (solid line), ao = 10 4 (dot- 
ted line), ao = 10 s (dashed line), in SI units. The vertical dashed 
lines denote the frequencies v\ , and V2 , where f 2 > v\ , which are 
used in section J3] 



material at r m i n becomes optically thin, so that the entire 
jet is optically thin to all higher frequencies. Hence, the ex- 
tent of the frequency range covered by the flat part of the 
spectrum for a given emission region can be determined by 
calculating the radii where t v — 1 as a function of frequency. 
This is described by the function, rth(v) oc v" 1 , and hence 
an emission region of size, r max /r m i n = 10 5 , should create 
a flat spectrum over the frequency range, ^ ma x/^min = 10 s . 
However, as figure [5] demonstrates, the flat region does not 
remain truly flat over the full span of this frequency range. 
This is because there is a smooth transition from the flat 
region of the spectrum to the optically thick and thin power- 
law tails. 

The observations of the GRS 1915+105 jet which are 
modelled in section ^| reveal a flux ratio between two fre- 
quencies that is close to, but not equal to, unity. If this 
slight deviation from a truly flat spectrum is interpreted as 
the effect of the optically thick cut-off slope, then for a fixed 
emission region size, the model parameter, ao, may be deter- 
mined by fitting the model to the observed flux ratio. The 
flux normalisation may then be used to determine Bo, and 
hence ko may also be extracted. 



either the electron density or the magnetic field strength 
throughout the entire length of the emission region. 

The length of the emission region is determined by the 
distance the ejected jet material has travelled in the time 
period defined by the flux oscillation. Emission from mate- 
rial that was ejected at earlier times, and hence is at greater 
radii, does not need to be modelled as it merely represents 
a constant flat background component which has no effect 
upon the time variability. Furthermore, the flux contribution 
of this background component to the flux peaks is insignif- 
icant, as this material will be located at radii where it is 
very optically thin to radiation of the observed frequency 
(see figure 0. 

From the numerical study of the previous section, the 
simplest mechanism to mathematically increase the flux nor- 
malisation of the flat spectrum would be a decrease in the 
magnetic field strength. However that requires ao to remain 
constant to ensure the flat spectrum still encompasses the 
same spectral range and hence the electron density would 
have to increase by an appropriate amount. It is not clear 
which physical process could cause such a fine-tuned change, 
and therefore the flux increase is best thought of as an in- 
crease in ao, due to an electron density increase and/or a 
magnetic field strength increase. 

In the absence of any clear physical expectations for 
time variability of the jet properties, and without loss of 
generality, we interpret the flux increase at the beginning 
of an oscillation as an increase in the electron density of 
the material injected into the jet. To simplify the model 
we begin and end with an injected electron density of zero, 
although in reality an increase from some quiescent electron 
density injection value is more likely. The resulting spectral 
behaviour of the model is the same, but in the latter case 
our results would just be superimposed upon a pre-existing 
steady flat spectrum. This underlying flat spectrum would 
be more extended then that which can vary in the observed 
time-scale, but that is not of any importance as it has not 
been observed. 

By allowing ko = ko (t) , we introduce a time dependence 
to the observed flux function. However, if the injected elec- 
tron density at the base of the jet varies with time, then it 
will also vary with radius, as there is a lag in the time for 
the injected electrons to reach greater radii along the jet. 
Therefore we introduce a transformed variable for t 



t'(t,r) = t- 



(13) 



3 TIME VARIABILITY IN FLAT 

SYNCHROTRON EMISSION SPECTRA 

In this section we investigate the possibility of explaining 
the time variability observed in jet sources with a time- 
dependent version of the BK79 model. The GRS 1915+105 
observations (see show that during an oscillation the 
ratio between the infrared and millimetre flux remains ap- 
proximately constant, indicating that the flat spectrum ap- 
parently persists throughout the full cycle of flux increase 
and subsequent decrease. Hence the time variability simply 
represents an oscillation in the flux level of the flat spec- 
trum, without any significant changes in its form over this 
frequency range. To achieve an increase in the flux normal- 
isation of the flat spectrum, in terms of the steady-state 
BK79 model, would require a near-simultaneous increase in 



Hence equation 1121 becomes 

_ 8.8 x TO" 18 1/2 / v n ■ ' 



with, 



/•i-maxW , N 3/2 , s 

J (^J [i-e-^)] drmJy, (14) 



(A/pc) 2 

f»-„ax(t) , v 3/2 



r v {t,r) = 1.5 x 10~ 12 k' (t') B 17/8 wo 

, - n -25/8 / v \ -25/8 

*y (gs) • ^ 

where r max (i) = rmin+Vjt. Here we have assumed a constant 
bulk velocity of the jet material, Vj. 

We implemented three variations of the electron density 
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Figure 3. Model light curves for V2 (solid line) and v\ (dashed line), where t>2 ^ > ^l? over a period of 2i max , for a) ao = 1.0 X 10 s and 
b) ao = 1.0 X 10 4 , with a constant value of injected electron density. 
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Figure 4. Model light curves for 1/2 (solid line) and v\ (dashed line), where v% ^> 1^1 } over a period of 2£ m axi 

for a) ao = 1.0 X 10 8 and 

b) ao = 1.0 X 10 4 , with a linear electron density injection profile. 




Figure 5. Model light curves for V2 (solid line) and v\ (dashed line), where V2 ^> z-'ij over a period of 2£ m ax> 

for a) ao = 1.0 X 10 8 and 

b) ao = 1.0 X 10 4 , with a Gaussian electron density injection profile. 



injection function, k' (t), with t expressed in units of t max - 
the time of the maximum injected electron density: 

1) A simple spontaneous increase in the injected elec- 
tron density, followed by a spontaneous decrease 



k' (t) = 



fc 

0.0 



< t < 1 
t > 1 



(16) 



This model will effectively produce the steady-state model 
spectra of figure|5|at, and only at, t = 1. 



2) A linear increase in the injected electron density, fol- 
lowed by a linear decrease 

*S(t) = *o(l-|t-l|). (17) 

3) A Gaussian electron density injection function 

k' (t) =fc e- 25( '- 1)2 . (18) 
Simulations were performed from t = to t = 2t max , 



where the values of t n 



and v\ are chosen such that the 
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emission region is the same length as the one used to create 
the steady-state spectrum in figure [5] The resulting light 
curves from two frequencies, j/2 and vi, where V2 > v\ (as 
defined in figure H, for each of the three injection functions 
are displayed in figures [3] [I] and [K] respectively. Each in- 
jection function light curve pair is shown for two different 
values of the model parameter ao, which is defined in section 

El 

As the increase in electron density propagates along the 
jet, a flat spectrum, which initially forms at the high fre- 
quency end, extends towards lower frequencies. The light 
curves peak at the time when the flat spectrum has ex- 
tended to include the frequency of observation. Hence, the 
light curve of the higher frequency emission always peaks 
before that of the lower frequency emission because it orig- 
inates from smaller radii within the jet. The profile of the 
light curve in figure H clearly demonstrates how a flat spec- 
trum is created as a function of time. The flat region of the 
spectrum is formed very rapidly at the high frequency end, 
but its progression towards lower frequencies slows down 
over time. When the electron injection is switched off, the 
high frequency end of the flat spectrum very rapidly disap- 
pears but the lower frequency region decays more gradually. 

For higher values of ao, the flat spectral region is shifted 
towards higher frequencies, so a given frequency peaks at a 
later time. Furthermore, for the lower frequencies, such as v\ 
in figure^ and for lower values of ao, the flat spectral region 
never reaches the frequency of observation, and so a slow flux 
increase, followed by a slow decay is observed. Conversely, 
as the observation frequency approaches the optically thin 
cut-off frequency, the observed light curve profile will tend 
towards the injection function profile. This may also be de- 
scribed by the optically thin radius for the frequency of ob- 
servation, v, approaching the base of the emission region, 
nh(f) — > r min . 

The linear and Gaussian injection functions both dis- 
play light curve profiles, shown in figures 2] and that tend 
towards the injection profile as r t h(^) — > ?"mm. The time lag 
observed between the flux peaks of the two frequencies can 
be reduced by decreasing ao, such that both frequencies lie 
closer to the, rapidly formed, high frequency end of the flat 
spectrum. This also has the effect of decreasing the peak flux 
ratio between the two frequencies, until both frequencies lie 
in the flat region where it becomes unity (assuming the emis- 
sion region is sufficiently large) . For sufficiently small values 
of ao the peak flux ratio will begin to decrease again, as the 
higher frequency emission becomes optically thin. 

For any fixed value of ao the predicted time lags, peak 
flux ratios, and flux normalisations are the same for both 
the linear and Gaussian injection functions, and this would 
be true for all injection functions with some gradual increase 
and decrease. Therefore for a given set of frequencies, and 
a given period of flux increase, t max , a value of ao may be 
determined from the observed time lags and peak flux ratios, 
independent of the exact functional form of k' (t). The flux 
normalisation then determines Bo, and hence ko may be 
extracted. This technique is employed in the case of the GRS 
1915+105 observations in ftOl 
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Figure 6. Infrared 2.2 fim data (crosses) from the IRCAM3 in- 
strument on UKIRT, and millimetre 1.3 mm data (diamonds) 
from the SCUBA instrument on JCMT. Taken from Fender & 
Pooley (2000). 



4 APPLICATION TO GRS 1915+105 

4.1 Observations 

GRS 1915+105 was observed simultaneously at two wave- 
lengths on 1999 May 20 (Fender & Pooley 2000); at 2.2 fim 
with the IRCAM3 instrument on the United Kingdom In- 
frared Telescope (UKIRT), and at 1.3 mm with the SCUBA 
instrument on the James Clerk Maxwell Telescope (JCMT). 
The resulting light curves recorded over a two hour period 
are shown, overlaid, in figure |S| The light curves show large 
amplitude quasi-periodic oscillations, each lasting ~ 1000 s, 
believed to be ejection events in the jet. Detailed views of 
the first and fourth peaks are shown in figure 

The infrared data from UKIRT has been dereddened 
with an assumed infrared A"-band extinction of Ak = 3.3 
mag. However, there is an uncertainty in the infrared flux 
values of at least 40 per cent (Fender et al. 1997) as the 
precise value for the absorption correction is not known. 
Furthermore we expect that the infrared flux will have some 
unknown background contribution due to the emission from 
regions of the GRS 1915+105 system other than the jet. The 
error values on individual data points are of the order of 5 
per cent (Fender, private communication). 

The millimetre data show flux peaks that occur almost 
simultaneously with the infrared peaks. The time lag be- 
tween the millimetre and infrared peaks is less than 50 s, 
though the sampling rate is not sufficient to provide a more 
exact determination. In the case of the second infrared peak 
there is no corresponding peak in the millimetre data. This 
peculiarity could prove useful in understanding the physics 
behind these oscillations (see H5.3I . 

4.2 Properties of GRS 1915+105 

We can use observations of GRS 1915+105 to set the fixed 
parameters of our model described in These parameters 
are summarised in tabled As explained in the inclination 
of the jet is implicitly assumed by the model to be i — 90° , 
which is in reasonable agreement with observed inclinations 
of i « 70° (Mirabel & Rodriguez 1994). The opening half- 
angle of the jet, #j, is observed to be less than 4° (Fender 
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Figure 7. Detailed view of a) 1st flux peak b) 4th flux peak in 
the data shown in figure 151 



Wo Dj Vj r min i max 

0.05 r 11 kpc 0.6 c 1.0 x 10 5 m 350 s 



Table 1. Fixed parameters of the jet model, determined by ob- 
servations of GRS 1915+105. 



et al. 1999), so we take wo ~ 0.05 ro. For the distance to 
GRS 19154-105, Dj, we assume 11 kpc (Fender et al. 1999). 
Following Kaiser et al. (2000) we adopt the value of Vj = 0.6c 
for the bulk velocity of the electron plasma in the jet. We 
neglect any relativistic effects on the emission spectra as 
these will be small for the adopted parameters (see H5.2H . 

We assume that the radius of the initial shock accelera- 
tion, r m i n , is where the relativistic electrons are first injected 
and thus is limited to the size of the last stable orbit for a 
non-rotating black hole, which is approximately three times 
the Schwarzschild radius, r s . Taking the lower mass limit 
for the GRS 1915+105 black hole of 10 M Q (Greiner, Cuby 
& McCaughrean 2001), gives the lowest plausible value for 
r m in = 3r s = 1.0 x 10 s m. The value of t maK is assumed to 
equal the rise time of the observed infrared flux peaks which 
is ~ 350 s. 




o L , , , i , , , , i , , , , i , , , , i , , , , i , , , , 

50 100 150 200 250 300 

Time-lag (s) 

Figure 8. The predicted time lag and flux ratios between the 
flux peaks at 1.4 X 10 14 Hz and at 2.3 X 10 11 Hz for 10 3 < a < 
3.5 X 10 s , where ro = 1.0 X 10 11 m. The grey rectangle shows 
the region where the predicted relationship would agree with the 
data. 

4.3 Light curve modelling 

The light curves seen from the GRS 1915+105 observations 
seem to be best fit by the Gaussian electron density injec- 
tion function. However, the light curve profile reveals a rise 
period that is apparently shorter than the decay period, sug- 
gesting a less symmetric functional form of k' (t) and hence 
disallowing a direct comparison between the model and data. 
Despite this, the model may still be used to study the time 
lags between the infrared and millimetre flux peaks, as this 
is a property of the light curve which is independent of the 
injection function profile (sec JJJ. 

The longer decay period may be suggestive of a delay 
between the time of maximum injected electron density and 
the time of the infrared flux peak. This may indicate the 
presence of a flat spectrum extending beyond the infrared 
into optical frequencies. However, for a fixed emission region 
size the extent of the flat spectrum beyond the infrared fre- 
quency of observation is solely determined by the value of ao, 
which is constrained by the observed time lags between flux 
peaks. Therefore, the light curve profiles are most probably 
due to an electron density injection profile effect. 

Using the model parameters of GRS 1915+105 shown in 
table Q we performed several simulations of the predicted 
light curves for the infrared (1.4 x 10 14 Hz) and the mil- 
limetre (2.3 x 10 11 Hz) fluxes with different values of the ao 
parameter. The resulting relationship between the time lag 
and the flux ratio between the flux peaks of the two frequen- 
cies is displayed in figure |H| Acceptable combinations of the 
time lag and the flux ratio allowed by the data are denoted 
by the grey rectangle. Obviously, the model does not agree 
with the data. Modelling smaller jet opening angles simply 
shifts the data points for a given value of ao down to a lower 
position upon the same relationship line. The flux ratio for 
a given time lag may be increased to agree with the data if 
the velocity of the jet was lowered significantly, but this is 
in strong disagreement with our understanding of the rela- 
tivistic jet. Finally, even substantial changes to the value of 
r m i n have only a small effect on the relationship. 

If our model is correct, then this finding suggests that 
the infrared flux data may not represent the true infrared 
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Figure 9. The model of the first infrared/millimetre peak as 
fitted to the observed peak fluxes and time lag, with the data 
overlaid using the same symbols as in figure [S] A better fit of 
the profile could be achieved by altering the electron density in- 
jection function, and by manipulating the balance between the 
background and extinction corrections. 



Figure 10. The predicted flux from 1.3 mm emission as a func- 
tion of the time lag between the infrared and millimetre flux 
peaks. Each line represents a different maximum injected elec- 
tron density, from the bottom line of k(r = r m i n ) = 10 16 (solid 
line), 10 22 , 10 28 , 10 34 , to 10 40 (dashed line) m~ 3 at the top. The 
time lags are solely determined by the values of ao- 



flux from the jet, due to non-jet related background emission 
and perhaps also an overestimate of the extinction (see also 
Fender et al. 1997). Therefore the true flux ratio may well 
be smaller. On the basis of this premise we now attempt to 
fit the model to the data, using the predicted flux ratio that 
corresponds to the observed time lag to determine the true 
infrared flux. We then check whether the thus determined 
infrared flux is within the observational uncertainty. If so, 
we can then determine the parameters fco and Bo from the 
model. 

4.4 Data fitting 

The observed time lag between the flux peaks of the two 
frequencies is approximately 25 s. The value of ao that pro- 
duces this time lag is 1.0 x 10 3 (with m = 1.0 x 10 11 m), 
for which our model predicts a flux ratio between the two 
peaks of 1.1. For the first infrared peak in the data, the cor- 
responding millimetre peak has a flux of ~ 340 mjy, so from 
our model we would then expect the infrared peak to have 
a flux of ~ 375 mjy. As the observed infrared peak has a 
flux of ~ 525 mjy in the data, we correct every data point 
by subtracting a hypothetical background emission compo- 
nent of 100 mjy, and then scale by a factor of 0.88 due to 
the uncertain extinction value. The same flux ratio may be 
achieved with a simple subtraction of a 150 mjy background 
component, however the quiescent infrared flux in-between 
flux peaks suggests that this value is too large. The reduced 
extinction is well within the 40 per cent flux uncertainty (see 

E3J- 

Figure |U] shows the model fit to the millimetre peak 
flux which obtained a value for the magnetic field strength 
at the radius of initial shock acceleration, r m i n = TO x 10 5 
m, of B(r — r m i n ) = 7.8 x 10~ 7 T, which corresponds, for 
the model value of ao, to a peak injected electron density of 
fc(r = r min ) = 5.4 x 10 40 m" 3 . 

This extreme value for the electron density cannot be 
avoided as is shown in figure [TU1 This plot shows the range 
of peak millimetre fluxes as a function of time lag for dif- 



ferent values of k(r = r m i n ), and B(r = r mln ). It is clearly 
impossible to obtain millimetre flux levels of the order of 
~ 100 mjy without extremely high electron densities, for 
time lags of less than 300 s. By using this time variability 
analysis we have avoided the need to rely upon the uncer- 
tain flux level of the infrared emission. Furthermore, this is 
not just a problem with this time dependent model, but is 
a fundamental problem of the BK79 model, as can be seen 
by applying the results of iJ2.2l to the observations of GRS 
1915+105. 

Finally, to maintain the emission of the highest fre- 
quency radiation at the outer edge of the jet, the elec- 
tron energy distribution must extend to a Lorentz factor 
of 7max — 10 7 for our fitted value of the magnetic field 
strength. Although this value is high there is some evidence, 
from AGN jet observations (Dermer & Atoyan 2002), that 
electrons in jets can be accelerated to such energies. How- 
ever, this high value of 7 max is a direct consequence of the 
low fitted value of the magnetic field strength, which would 
be higher in a model that does not require extremely high 
electron densities. 



5 DISCUSSION 
5.1 Summary 

In this paper we have demonstrated how the BK79 model 
produces a flat synchrotron spectrum from a conical jet ge- 
ometry following the assumption that the emitting electron 
plasma is free of energy losses as it expands away from its ori- 
gin. We have shown that the frequency range covered by this 
flat spectrum depends upon the radial extent of the emission 
region with respect to the radius where the electron plasma 
becomes optically thin to radiation at each frequency. This 
radius is determined by the value and radial dependence of 
the magnetic field strength and electron density within the 
jet. From observations of the frequency range that the flat 
spectrum extends over, and the observed flux from these fre- 
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quencies, both the magnetic field strength and the electron 
density in the jet may be individually determined. 

We also adapted the BK79 model to include time vari- 
ability. This not only allows comparisons with light curve ob- 
servations, and places stronger demands upon the hitherto 
unexplained re-energisation process, but also places more 
stringent constraints upon the model. Prior to the inclu- 
sion of time variability the extent of the emission region 
had to be assumed. The time scale of variability now sets 
the spatial extent of the emission region. Furthermore the 
model parameters may be derived from observation of the 
time lag between peaks at each frequency rather than the 
less reliable measure of the flux ratio between the peaks at 
each frequency. Time lags are caused by an optical depth 
effect that occurs whether or not the flat spectrum extends 
beyond the observed frequency range, and hence time lags 
are a more reliable indicator of the physics of the system 
than the flux ratio. The predicted time lag would be the 
same even with an inverted spectrum that fully covers the 
observed frequency range. 

This time variability model was then applied to obser- 
vations of the GRS 1915+105 microquasar, which revealed 
a discrepancy between the model light curve profile and 
the observed profile. Furthermore, the model calculated for 
the known parameters of GRS 1915+105 required extremely 
large electron densities to explain the observed flux. As we 
believe that the system parameters are sufficiently well con- 
strained, the only conclusion must be that the model itself 
is flawed. 

5.2 On the flux discrepancy 

As demonstrated in figure fTOl we have shown that the ob- 
served millimetre flux levels are only predicted by the mod- 
els that also predict time lags, between the infrared and 
millimetre peaks, which are many times greater than those 
observed (for physically justifiable electron density values). 
We believe that the observed millimetre flux levels are an 
accurate representation of the true millimetre jet emission, 
with no significant background contamination, as the emis- 
sion drops to zero between flux peaks. Therefore, the model 
is fundamentally incapable of predicting millimetre flux lev- 
els of the order of 100 mjy, together with time lags between 
flux peaks that are less than 100 seconds, which is what has 
indisputably been observed. 

The bulk gas flow within the jet is relativistic. There- 
fore, if relativistic effects can cause the observed time lag to 
be considerably shorter than the time lag in the rest frame 
of the jet material, then the observed flux will be predicted 
from a lower electron density. Although time dilation has the 
opposite effect, relativistic Doppler shifts can in principle 
cause this effect. However, in the case of the GRS 1915+105 
jet, the inclination of the jet axis to our line-of-sight is well 
constrained to ~ 70°, and for such inclinations the max- 
imum effect on the observed flux from relativistic Doppler 
shift and Doppler beaming is only 20 per cent. To explain the 
observations with believable electron densities the predicted 
flux needs to be 10 4 times greater. Such a difference cannot 
easily be rectified by geometrical effects, nor by the current 
uncertainty in the distance to GRS 1915+105. Furthermore, 
the simplification in using the equivalent line-of-sight opti- 
cal path length for a jet at an inclination of 90° does not 



significantly effect the results. Only for substantially smaller 
viewing angles and highly relativistic jet velocities do rela- 
tivistic effects lead to a significant shortening of observed 
time lags compared to the jet rest frame. 

The observations are almost certainly of synchrotron 
emission, as thermal emission would require even greater 
electron densities (Dhawan et al. 2000) , and polarisation ob- 
servations agree with this assertion (e.g. Fender et al. 2002). 
The flux discrepancy cannot be explained by the presence 
of the more extended jet beyond the region of the time vari- 
ability, as the time variability still has to explain a large flux 
increase. 

In conclusion it is impossible to explain synchrotron 
emission of this strength from a BK79 type jet with justifi- 
able electron densities, when restricted to the observed time 
lags. 

5.3 Is the flat spectrum flat? 

From a sample of just two frequencies we cannot conclusively 
say that the observed flux ratio is due to a flat spectrum 
which does not completely extend to the lower frequency. 
The other possibility is that the 'flat' region of the spectrum 
extends to cover both of the observed frequencies and has a 
slightly positive slope. Such a scenario would naturally arise 
if the energy losses in the jet due to adiabatic decompres- 
sion were partially included. However, if adiabatic losses are 
fully included then the resulting inverted spectrum would 
have a slope that is significantly steeper than that observed. 
The steady-state jet model of Falcke (1996) includes partial 
energy losses within the jet, and, for the GRS 1915+105 jet 
inclination angle, predicts a spectral slope of a m 0.2. The 
observations allow a spectral slope of 0.0 < a < 0.12, due to 
the uncertainty on the infrared flux, so we cannot exclude 
the possibility of a fully extended, truly flat spectrum. 

Hence, there are two possible scenarios to explain a non- 
unity flux ratio between the flux peaks of the two observed 
frequencies. Either the lower frequency observation is of the 
optically thick region (or the transition to this region) of the 
spectrum from a BK79 type jet, or it is of the 'flat' region 
of a jet with energy losses. There exists a simple observable 
difference between these two scenarios; whereas for the for- 
mer case a relationship exists between the time lag and the 
flux ratio, for the latter case these two parameters are inde- 
pendent of each other. Although the loss of the relationship 
between the time lag and flux ratio will not in itself affect 
the predicted electron densities, its implication of energy 
losses within the jet could lead to a slightly lower predicted 
electron density. 

The missing millimetre flux peak to coincide with the 
second observed infrared peak does seem to suggest that for 
the second ejection event the flat spectrum does not extend 
all the way down to the millimetre waveband. However, this 
explanation is applicable to either of the scenarios discussed 
above, and is perhaps more likely than a shift in the spectral 
slope. 

5.4 Conclusion 

We have developed a time-dependent version of the par- 
tially self-absorbed jet model of Blandford & Konigl (1979). 
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As expected from the original steady-state model, the time- 
dependent model gives rise to flat broadband spectra but 
cannot explain the large flux variations observed in GRS 
1915+105 at millimetre and infrared frequencies without 
invoking unrealistically high electron densities. This result 
also holds true if the time variability is ignored. In this case 
the parameters of the steady-state model are determined 
through the observed flux ratio of the two frequencies at 
any given time, rather than from the observed time lag be- 
tween flux peaks of the two frequencies. Therefore either this 
variable emission does not originate from the jet of GRS 
1915+105, or an alternative model is required to explain 
microquasar jet emission. 

Observations of the jet spectrum and time variability 
at different frequencies will help resolve the issue of whether 
the spectrum is flat, and whether we are observing optical 
depth effects. Observing a relationship between the flux ratio 
and the time lag would confirm that the flux ratio is due to 
an optical depth effect rather than due to a mildly inverted 
spectrum. 

Recent observations of the infrared flaring behaviour 
of GRS 1915+105 (Rothstein & Eikenberry 2002) have il- 
lustrated the complex nature of the time variability in the 
emission from relativistic jets. Hence, more complex models 
will ultimately be required to fully explain the observations. 
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